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Abstract 

Using the HI detector at HERA, charged particle muhiphcity distributions in deep in- 
elastic e^p scattering have been measured over a large kinematical region. The evolution 
with W and of the multiplicity distribution and of the multiplicity moments in pseu- 
dorapidity domains of varying size is studied in the current fragmentation region of the 
hadronic centre-of-mass frame. The results are compared with data from fixed target 
lepton-nucleon interactions, e^e~ annihilations and hadron-hadron collisions as well as 
with expectations from QCD based parton models. Fits to the Negative Binomial and 
Lognormal distributions are presented. 
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1 Introduction 



The multiplicity distribution of hadrons produced in high energy interactions is one of the 
basic measures characterising multiparticle final states. The fluctuation pattern of the number 
of particles produced in a given domain of phase space reveals the nature of the correlations 
among the hadrons and is, therefore, sensitive to the dynamics of the process. The multiplicity 
distribution of charged hadrons has been measured for the full variety of collision processes 
from e~^e~ annihilations to nucleus-nucleus collisions and the available data cover a wide energy 
range Numerous theoretical studies have been devoted to the subject [§, starting 

with early investigations by Heisenberg and Fermi 0. 

Whereas the total event multiplicity remains of considerable importance, interest has shifted 
with time towards studies of the multiplicity distribution in subdomains of phase space. In these 
restricted domains, global conservation constraints are minimised and dynamical correlation 
effects better revealed 0,11. 

In this paper we present results on the multiplicity distribution of charged hadrons pro- 
duced in the current fragmentation region of deep inelastic scattering (DIS) e~^p collisions. The 
analysis is based on data accumulated by the HI detector at HERA in 1994, corresponding to 
an integrated luminosity of 1.3 pb^^. 

The multiplicity distribution and its moments, fully corrected for detector effects, are stud- 
ied in the virtual-boson proton {'y*p ) rest system, the hadronic centre-of-mass frame. The 
measurements are obtained in subdomains of pseudorapidity spaceQ, both as a function of W, 
the total hadronic centre-of-mass energy, and Q^, the negative four- momentum transfer squared 
in the DIS process. Results from deep inelastic e~^p interactions from the HERA collider on 
the charged multiplicity distribution and on the mean multiplicity in the current region of 
the Breit-frame |]T0[ have already been published. 

Here, HI results are compared with those obtained in DIS fixed target experiments at much 
lower energy, with hadroproduction in e~^e~ annihilation reactions, with data from hadron- 
hadron collisions as well as with the expectations of QCD based parton shower models and 
other more phenomenological approaches. 



2 Definitions and phenomenology 
2.1 Correlations 

The set of probabilities P„ for various numbers of charged hadrons {n) to be produced in a 
given region of phase space is known as the multiplicity distribution. The continued interest 
in multiplicity distributions rests on the observation that fluctuations in "counting statistics" 
are a direct measure of the strength of correlations among the objects being counted. In high 
energy physics, this was first exploited by Mueller [|lT| in the formulation of the concept of short 



range order. Since then, correlations have been extensively used as probes of the interaction 
and hadronisation dynamics 0]. 

The connection between the multiplicity distribution and particle correlations is made ex- 
plicit in the relation between the factorial moment of order q, Rg, of the multiplicity distribu- 
tion in a phase space domain (say, rapidity A) and the g-particle inclusive momentum density 

^Pseudorapidity is defined as rj* = — lntan|, with 9 the angle between the hadron momentum and the 
direction of the virtual photon in the rest system. The current hemisphere is defined as the region of 
positive 77*. 
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Pq(l, 2, . . . , g) which for identical particles takes the form |TT|: 



Rq = ^n{n-l)...{n-q + l)Pn= ... Pg{l,2, . . . ,q) dyi . . . dyq. (1) 

Q J A J A 

The functions Pg contain, besides genuine dynamical or kinematical correlations, contributions 
from "random coincidences" in the region A. The latter are eliminated by considering the 
("connected") correlation functions 2, . . . , g), familiar from statistical physics. Integrated 
over a domain of size A they define the factorial cumulants (or Mueller moments) Kg of the 
multiplicity distribution For the normalised quantities Rg = Rg/ {nY and Kg = Kg/ {nY 
the following relations hold: 

K2 = R2-I, (2) 
Ks = R3-3R2 + 2; (3) 

where (n) is the mean multiplicity in A. 
2.2 Scaling 

Koba, Nielsen and Olesen (KNO) ||T^ have studied the limit of the multiplicity distribution in 
full phase space for n ^ 00 and {n) 00 with z = n/ (n) fixed. In this case one obtains the 
KNO form 

(n)P„~^(;2). (4) 

The function "^{z) was shown to become asymptotically independent of the total energy if 
Feynman scaling [^] is satisfied. A proper mathematical reformulation of multiplicity scaling 
for discrete distributions, valid also at finite energies, was later given by Golokhvastov and 
is known as KNO-G scaling. 

Exact KNO scaling implies that, besides "^{z), the factorial moments Rg as well as the 
moments Cg = {n'^) / {nY are energy independent. 



Scaling, either in KNO or in KNO-G form, is experimentally well-established [[L5ip!6| for 
the full phase space and single hemisphere multiplicity distributions in e~^e~ annihilations, in 
DIS lepton-hadron interactions and in hadron-hadron collisions [|l7l, except at the highest SPS 



collider energies [|T8],0. This is remarkable since Feynman scaling is strongly violated in all 



these processes. However, already in 1970, Polyakov |2y] derived the KNO scaling law for 
e'^e~ hadronic final states within a broad class of field theories. This pre-QCD model is based 
on a conformal invariance principle which can be reformulated in terms of a scale-invariant 
stochastic branching process. An energy independent coupling constant is assumed and the 
mean multiplicity rises as a power of the energy. The Polyakov derivation, and subsequent 
work |^0|, demonstrates that KNO scaling should hold at least approximately in any model 
based on similar principles. 

2.3 Parametric models 

A much studied distribution is the Negative Binomial Distribution (NBD) defined as 
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with parameters k (or 1/k) and n. The average (n) and the dispersion D of the NBD are 
related to the two parameters by 



(n) 



n 



1 1 
- + 



n 



k 



(6) 



For 1/k (H) reduces to a Poissonian. n. The parameter is equal to the integrated 
second order correlation function K2 in the studied phase space domain. 

Many phenomenological models for hadropro duct ion predict multiplicity distributions of 
a Negative Binomial form. They are reviewed in 1^,^. In QCD, the NBD is obtained as a 



solution in Leading-Log Approximation for the gluon multiplicity distribution in a quark jet 



Within the framework of the Modified Leading-Log Approximation (MLLA) and Local Parton- 
Hadron Duality (LPHD) ||2^ , the lower order factorial moments behave approximately as those 
of the NBD ||,||,||]. 

It is well-established experimentally in a large variety of collision processes that multiplicity 
distributions, in full phase space as well as in restricted phase space domains, are approxi- 
mately of Negative Binomial form. However, deviations are observed in high statistics e~^e~ 
experiments for final states with hard jets ||27|,|28[| and in hadron-hadron interactions |jT9|. In 



fixed target DIS lepton-hadron and lepton-nucleus collisions the NBD adequately describes the 
multiplicity distributions in the full phase space, in limited (pseudo-) rapidity intervals, as well 
as in the full current and target hemispheres [@,|2S|. 

A multiplicity distribution exhibiting explicit KNO-G scaling has recently been derived by 
assuming a scale-invariant bi-variate multiplicative branching mechanism as the basis of mul- 
tihadron production [|16[. Application of the central limit theorem leads to a scaling function 
of Lognormal form. In this model the multiplicity distribution P„ is related to a continuous 
probability density f{n) and defined as P„ = J^~^^ f{n)dh, where f{h) is the Lognormal distri- 
bution (LND). The mean continuous multiplicity (n) is approximately given by (n) = (n) + 0.5. 
Following W^ one finds 



("+!)/("> 
n/{h) 



N 



2TTa z + c 



exp 



[ln{z + c) - fif 
2^ 



dz. 



(7) 



The integrand in defines the Lognormal scaling function in KNO-G form. It depends on 
two parameters. Here, z = h/ (n) is the scaled continuous multiplicity; A^, a, /i and c are 
parameters of which only two are independent due to normalisation conditions. In fits to data, 
correlations between the parameters are reduced if d and c are used as free parameters. In 
terms of o"^ and fi these are given by 



\ 



In 



d 



1 + 



+ 1 



and fx = ln(c + 1) — cr^/2. 



(8) 



The parameter d is equal to the dispersion of the scaling function. The expression for 
Pn also depends on (n) or (n). The latter fit parameter is denoted by m in Table |^. Exact 
scaling implies that d and c are energy independent. In addition, the generalised dispersions 

= [((^ " i''^))'^)]^^'' (-0 = ^2) satisfy a generalised Wroblewski relation Finally, the 



mean continuous multiplicity is found to grow as a power of the energy. Comparisons of the 
Lognormal distribution with data from e~^e~ annihilation, v{V)p and pp collisions can be found 
in |8|,|1|,|1],||,|3[. 
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3 Experimental procedure 



3.1 The experiment 

The experiment was carried out with the HI detector at the HERA storage ring at DESY. 
The data were collected during the 1994 running period when positrons, with an incident energy 
of 27.5 GeV, collided with protons with an energy of 820 GeV. The following briefly describes 
the detector components most relevant to this analysis. 

The energy of the scattered positrons is measured with a liquid argon (LAr) calorimeter and 
a "backward" electromagnetic lead-scintillator calorimeter (BEMC). The LAr calorimeter P3 



extends over the polar angular range A° < 9 < 153° with full azimuthal coverage, where 
6 is defined with respect to the proton beam direction (positive z axis). It consists of an 
electromagnetic section with lead absorbers and a hadronic section with steel absorbers. 

The BEMC covers the polar angular range 151° < 6 < 177°. A principal task of the BEMC 
is to trigger on and measure scattered positrons in DIS processes with values ranging from 
5 to 100 GeVl The BEMC energy scale for positrons is known to an accuracy of 1%. 

The central tracking detectors are a hybrid of inner and outer cylindrical jet chambers (C JCl 
and CJC2), 2;-drift chambers and proportional chambers. The latter provide a fast signal and 
allow HI to trigger on tracks which originate from the z range expected for e~^p collisions. 

The jet chambers, mounted concentrically around the beam line, provide particle charge and 
momentum measurements from track curvature and cover the angular interval 15° < 6 < 165°. 
Up to 56 space points can be measured for non-curling tracks. 

The calorimeters and central trackers are surrounded by a superconducting solenoid pro- 
viding a uniform magnetic field of 1.15 T parallel to the beam axis in the tracking region. 

A backward proportional chamber (BPC), situated immediately in front of the BEMC and 
with an angular acceptance of 155.5° < 9 < 174.5°, serves to measure the impact point of 
the scattered positron and to confirm that the particle entering the BEMC is charged. Using 
information from the BPC, the BEMC and the reconstructed event vertex, the polar angle of 
the scattered positron can be determined to better than 1 mrad. 

Behind the BEMC is a time-of-fiight system with a time resolution of about 1 ns. This 
enables rejection of background events from interactions of protons in the material upstream 
of the HI detector. 

Positrons and photons emitted at very small angles to the positron beam direction are 
detected in two electromagnetic calorimeters located at 33 m and 103 m from the interaction 
point in the positron direction. Designed to measure luminosity by detecting 67 coincidences 
from quasi-elastic radiative e+p collisions, they are in addition used for studies of background 
arising from photoproduction. 



3.2 Event and track selection 

In this analysis two event samples are considered. The "low-Q^" sample consists of events in 
which the scattered positron is detected in the BEMC and is limited to a range in from 10 to 
80 GeV^. The events of the "high-Q^" sample have a positron detected in the LAr calorimeter 
and have > 200 GeV^. The event kinematics are reconstructed using the "lepton-only" 
method based on information from the reconstructed positron and the event vertex. 

Neutral current DIS events are selected by demanding a well-reconstructed scattered posi- 
tron with an energy larger than 12 GeV. This ensures that the remaining photoproduction 
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background comprises less than 1% of the DIS data sample. 

Various other sources of background are further reduced by appropriate selections. To 
exclude events with large QED radiative effects and to ensure substantial hadronic energy flow 
in the detector, the invariant mass squared, W"^, of the hadronic system, determined from 
energy clusters in the calorimeter, is required to be larger than 3000 GeV^. 

For comparisons with e~^e~ annihilation and non-single-diffractive data in hadron-hadron 



collisions, so-called "rapidity-gap" events are removed from the event sample. To that end, 
events for which the energy deposited within the polar angular range 4.4° < 9 < 15° is lower 
than 0.5 GeV are excluded. The same selection is applied to the various Monte Carlo generated 
event samples used in the analysis. 

An event vertex, reconstructed from tracks in the Central Tracker and located within ±30 cm 
of the mean vertex z position, is required to reject beam-induced background and to permit a re- 
liable determination of the kinematic variables. Remaining background is rejected by requiring 
no veto from the time of flight system. 

For this analysis, only charged particle tracks in the Central Tracker are used, after rejection 
of a possible positron track. Tracks are required to originate from the primary vertex and to 
lie within the polar angular range 15° < 6 < 155°. The higher limit eliminates the scattered 
positron in the low-Q^ sample. In the high-Q^ data the positron track is identified with a 
track-cluster linking algorithm. 

With the above selection the overall efficiency for finding a genuine primary track is better 
than 95%. This estimate is based on Monte Carlo simulations and cross-checked in a visual 
scan of real and simulated events. From this scan we further conclude that the tracker efficiency 
is simulated to an accuracy of better than 2%. 

The contamination of genuine primary tracks by tracks that have been fitted to the primary 
vertex, but originate from decay products of short lived particles, from secondary interactions 
and from photon conversions, is estimated to be 10%, < 1% and ~ 1%, respectively. No 
spurious tracks due to noise hits have been found. 

To investigate the sensitivity of our results to the track selection the analysis has been 
repeated after changing various track quality criteria, such as the number of hits per track, 
the cuts on minimal radial track length, the position in the CJC of the first measured point of 
a track, the minimum laboratory transverse momentum and the maximum distance of closest 
approach to the event vertex in the transverse plane. The results reported below are found 
to be stable against all these variations. Residual differences are included in the systematic 



uncertainties, further discussed in Sect. 3.5 



Within the considered kinematical regions, 53109 low-Q^ and 1576 high-Q^ events satisfied 
the above selection criteria, corresponding to a total integrated luminosity of 1.3 pb~^. The 
characteristics of the various event samples retained for further analysis are detailed in Table |l[ 

3.3 Models for hadronic final states 

The hadronic final state of neutral current DIS events is modelled using several Monte Carlo 
event generators. 

The event sample labelled as "DJANGO 6.0" has been generated with the DJANGO pro- 
gram It is based on HERACLES |]38| for the electroweak interaction and on the LEPTO 



program to model the hadronic final state. HERACLES includes first order radiative QED 



corrections and the simulation of real Bremsstrahlung photons. The structure functions used 



are based on the GRV parameterisation |^ and describe the HERA F2 results well. LEPTO 



8 



is used with the colour dipole model as implemented in ARIADNE pT| to model QCD par 



ton cascades. The hadronisation phase is modelled with JETSET and is based on string 
fragmentation . 



The sample hereafter labelled as "MEPS 6.4" has also been generated with LEPTO but dif- 
fers in the treatment of the parton cascade by using first order matrix elements and leading-log 
parton showers. It does not include initial or final state QED radiation. For the determination 
of systematic errors connected with the topology of the hadronic final state earlier versions of 
the MEPS generator as well as HERWIG Q have been used. 

For the events generated with the models listed, the detector response was simulated using 



a program based on the GEANT |^ package. The detector simulation has been checked by 
comparing in detail distributions of track quality estimators with the data. 

The DJANGO 6.0 sample, comprising nearly 120K events after all selections, is used here 
for data correction and to unfold the raw multiplicity distributions. The other Monte Carlo 
generated samples serve to cross-check the multiplicity unfolding procedure and to determine 
systematic uncertainties. 

3.4 Data correction 

The raw multiplicity distribution in a given region of (W, Q"^) and pseudorapidity needs to be 
corrected for several effects. These include loss of events and particles due to limited geometrical 
acceptance and resolution of the tracking system, limited track finding efficiency, contamination 
by tracks from particle decays and interactions in the material of the detector which are assigned 
to the primary vertex, and also QED initial state radiation which affects the event kinematics. 
The results presented in Sect. 4 are corrected for all the above-mentioned effects using the 
DJANGO 6.0 event sample. 

Correction factors are obtained from Monte Carlo simulation by comparing the "true" gen- 
erated distributions before the detector simulation with the "observed" distributions after this 
simulation followed by the same reconstruction, selection and analysis as the real data. The 
"true" distributions do not include the charged decay products of K^, A, A and from weakly 
decaying particles with lifetime larger than 8 ■ 10~^ s. 



The multiplicity unfolding method used here is similar to that first employed by TASSO [46 



and subsequently adopted in various other analyses ||,|2^,|3T],^. Let T^^^^{m) be the "true" 



number of events produced by the Monte Carlo generator with "true" primary charged hadron 
multiplicity m, and i?^^(m) the number of events with "true" multiplicity m observed after 
full detector simulation, kinematical reconstruction and event selection. The number of events 
with m generated tracks and n tracks accepted after simulation and track selection, divided by 
the total number of events observed with n accepted tracks, R^^^\n), determines a response 
matrix with elements A{m,n). The element A{m,n) is the fraction of observed events with n 
accepted tracks that had "true" multiplicity m. The observed multiplicity distribution Rf^^^{n) 
and the "true" multiplicity distribution T^^^^{m) are then related by 

TtTeim) = E n) Rt':{n). (9) 

The first factor at the right-hand side of corrects for any complete loss of events caused e.g. 
by inefficiencies and migration outside the investigated kinematical region due to inaccurate 
reconstruction of the event kinematics. 
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One of the underlying assumptions of the method is that the relative frequencies occurring in 
the response matrix derived from the Monte Carlo generator are the same as in the real data. 
Otherwise, the unfolded multiplicity distribution will be biased towards the generator input 
distribution. To reduce this bias it has been found necessary to use an iterative procedure 
whereby the predicted multiplicity distribution at the generator level is reweighted using a 
previous approximation of the unfolded data multiplicity, until convergence is reached^. This 
is illustrated in Fig. |l| which shows the raw multiplicity distributions in the interval 1 < rj* < 
5 for intervals in in a linear (upper figures) and logarithmic scale (lower figures). The 
dotted histograms are results from the DJANGO program after detector simulation and deviate 
appreciably from the real data. The result of the iterative reweighting method is shown as the 
solid histogram and reproduces the data very well. 

As a cross-check of the measurement of the means and dispersions of the distributions ob- 
tained by the matrix method, and to allow the study of certain multiplicity distributions given 
in parametric form, a second method is used. Instead of correcting the observed multiplic- 
ity distribution, a theoretical distribution T/^'j^g(m, a), depending on the parameter set a, is 
transformed into a raw reconstructed distribution taking all detector effects into account. The 
function Tl^^^{m,a) represents the true multiplicity distribution for ideal non-radiative deep 
inelastic collisions. The transformed distribution i?*g^(n) is then compared to the measured 
one and the parameters a determined with a standard minimum procedure. 

A more detailed discussion of the merits and limitations of this method is given in p9|,^ 



In the present paper it is used to make comparisons with the Negative Binomial and Lognormal 
distributions. 



3.5 Systematic errors 

In this section we discuss various sources of possible systematic uncertainties. Each possible 
systematic effect is independently varied and the analysis repeated as outlined in the previous 
section. The difference in the final result, relative to the quoted result based on DJANGO 6.0, 
is taken as one contribution to the systematic uncertainty. The errors from all sources are 
added in quadrature. 

For illustration, the change in the mean charged multiplicity due to various effects are 
summarised in Table ^ Results are given for both the full current hemisphere {vj* > 0) where 
the corrections are largest and the pseudorapidity range 1 < r^* < 3, where the detector 
acceptance is high. 

• Photoproduction events that survive the DIS selection will have a fake positron detected in 
the BEMC. Above a reconstructed scattered lepton energy E'^ of 14 GeV the contribution 
of photoproduction events is practically zero. For 12 < E[ < lA GeV the background 
contribution is estimated to be 5% [^. This means that photoproduction events only 
contribute in the highest W interval. To exclude events in this interval the E'^ cut is raised 
from 12 to 14 GeV, a smaller radial extension of the BEMC energy cluster is demanded 
and events with a positron detected in the electron tagger are rejected. This results, 
however, in a negligible change of the mean multiplicity in the highest W interval. 

• The BEMC energy scale is known to an accuracy of 1%. Decreasing the energy scale by 
its error reduces the mean multiplicity by a maximum of ~ 3%. 

^Convergence was reached in all kinematical regions after 3 iterations, at the most. 
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• The track selection criteria ensure a high efficiency (above 95% inside the acceptance 
region). To restrict tracks to a region of even higher efficiency additional kinematical 
selections are introduced. The polar angle ^ of a track is required to be larger than 22° 
and the transverse momentum px to be larger than 0.15 GeV/c. After correction, this 
decreases the mean multiplicity in the full current hemisphere by about 2%. However, 
exclusion of tracks from the forward Central Tracker region or with low px enhances the 
sensitivity to the type of Monte Carlo generator used. The contribution to the total 
systematic error is, therefore, estimated by applying additional selections on the track 
quality, requiring more than 10 hits, a track length of more than 10 cm and a start point 
in CJCl. With these selections, the fraction of tracks for which track segments in the 
inner and outer CJC are not linked, and therefore counted double, is kept below 0.3% 
and reproduced by the Monte Carlo simulation within 30%. 

• The multiplicity unfolding is performed iteratively to reduce the dependence on the gen- 
erated Monte Carlo multiplicity distribution. To estimate the importance of any residual 
systematic bias on the mean and dispersion of the multiplicity distribution, a Negative 
Binomial distribution, smeared for detector effects, is fitted to the reconstructed data 
multiplicity distribution. The effect on the determination of the mean multiplicity is well 
below 2%. In addition, {n) has been estimated from the fully corrected single-inclusive 
charged particle rj* spectrum in each W interval, using the standard bin-by-bin correction 
method. The (generally small) difference between extremes is taken as one contribution 
to the overall systematic error. 

• The predictions for (n) of the Monte Carlo generators used here differ by 17%, at most, in 
the pseudorapidity domain < 77* < 1. This region lies partly outside the Central Tracker 
acceptance at low W . Since a change in generator multiplicity outside the acceptance 
region does not affect the reconstructed distribution, a generator dependent bias cannot 
be removed by reweighting the input distribution. A systematic uncertainty has been 
assigned to the results pertaining to the full current hemisphere {rj* > 0) by taking the 
maximum difference between results derived from different event generators. Table 2 
illustrates the maximum size of the effect obtained from the generator LEPTO 6.1. 

The contribution from decay products of Kg, A, A etc. is subtracted via the unfolding 
procedure. However, recent studies at HERA indicate that the Ka production rate may be 



overestimated by about 10% in the models used for correction. Such a difference leads 

to a systematic underestimation of the mean multiplicity of at most 0.4%. 

From event simulation it is estimated that about 0.1 tracks per event assigned to the pri- 
mary vertex are, in fact, due to photon conversions in the detector material. To account for 
possible differences between the event simulation and the real detector response an additional 
1% systematic error on the mean multiplicity is assumed. 



4 Results 

The multiplicity distribution has been measured, in the kinematic regions in W and Ugted m 
Table 1, for charged particles with pseudorapidity in the domains I < r]* < rj* with rj* = 2, 3, 4, 5 
and in intervals of unit pseudorapidity centered at rj* = 2.5, 3.5, 4.5, as well as for the full current 
hemisphere defined as the domain 77* > 0. The data are integrated over particle centre-of-mass 
transverse momentum to allow comparison with other experiments. 
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The corrected multiplicity distributions, in four W intervals and four r]* domains, are given 
in Table 3. For each multiplicity, n, the relative frequency, P„, the statistical uncertainty 
(first) and the systematic uncertainty (second) are given. The statistical uncertainties take 
into account the finite number of events in the data and in the generated samples. They 
are calculated by means of a Monte Carlo sampling procedure which ensures propagation of 
sampling fluctuations and statistical correlations into the statistical error of the measurement. 
The same technique is employed to determine statistical errors on moments of the multiplicity 
distribution and other derived quantities. It must be noted that the P„ values given in the 
table are correlated for nearby values of n due to the method of correction. 

Table 4 summarises corrected results in various pseudorapidity domains for the mean charged 
multiplicity (n), the generalised dispersions Dq, the normalised multiplicity moments Cg {q = 
2, 3, 4) and the factorial moments R2, -R3. The table also lists values for the third order factorial 
cumulant for those intervals where a statistically significant measurement is possible. 

Estimates of the first two moments (i.e. (n) and D) of each multiplicity distribution have 
also been obtained using the parametric method, described in Sect. 3.4, based on fits of the 
Negative Binomial and Lognormal distributions. The results from this method and from the 
matrix unfolding are found to be consistent. 

Not all of the statistical moments listed in Table 4 will be discussed in detail. They have 
been measured in many experiments and are included for reference purposes. 



4.1 Current hemisphere multipUcity moments 

To allow comparison with lepton-nucleon fixed target data and single hemisphere data from 
e~^e~ experiments, we present in this section results on the moments of the charged particle 
multiplicity distribution in the full current hemisphere t]* > 0. The limited experimental 
acceptance in the interval < 77* < 1, in particular for the range 80 < 1^ < 115 GeV, renders 
the HI data in this region more sensitive to the Monte Carlo generators used in the correction 
procedure than in other pseudorapidity domains. The quoted systematic errors reflect this 
additional uncertainty^. 

Data on multiplicity distributions in a single hemisphere and in restricted domains of rapid- 
ity are available from several e~^e~ experiments p7| , |31| , ^,p2|. To compare DIS data with that 
for e~^e~ annihilations we have chosen to use the e~^e~ JETSET parton shower model, with pa- 



rameter settings as used by the DELPHI collaboration ||53[ . This model reproduces in detail the 
multiplicity distributions in the full phase space and in restricted intervals of rapidity [p7| ,p8|. 

The JETSET e~^e~ predictions presented below (labelled as "JETSET e~^e~" in the figures) 
are obtained for a mixture of "primary" light quark pairs only. Decays of Kg, A and A are 
treated as in the DIS data and in the simulations, thereby avoiding experiment and energy 
dependent corrections to published dataQ. The contribution to (n) from charmed quark frag- 
mentation, estimated at LEP to be a factor 1.02 ± 0.03 larger than for light quarks is 
neglected. The model predictions are also used in domains of pseudorapidity where no direct 
e~^e~ measurements exist. For these, the conclusions should be treated with caution. 



DifFractive events are removed from the HI DIS sample by rejecting events with a "rapidity gap" . Not 
removing these events would result in a 3% decrease of the mean multiplicity. This effect is not included in the 
quoted systematical errors. 

"^We have verified that the JETSET e+e^ predictions, with five active flavours, reproduce the published data 
in the PETRA-LEP energy range, including recent measurements at 130 GeV fs^ . 
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4.1.1 Mean Charged Multiplicity 

Figure 2a shows the mean charged multiplicity for rj* > 0, measured by this and other lepton- 
nucleon experiments |^,^|5^,0 as a function of W. 



Here and in the following, unless stated otherwise, the total errors are the overall uncer- 
tainties computed by adding the statistical and systematic errors in quadrature. When two 
error bars are displayed the inner error bar is the statistical error and the outer one is the total 
error. For data from other experiments we have used the published systematic errors, whenever 
available. Otherwise a systematic error of 5% is assumed. 

In the W range covered by HERA, (n) is compatible with a linear increase with InVT. 
Combined with the data at lower energy^ it is evident, however, that the mean multiplicity 
increases faster than In W. The HERA data confirm, for the first time in DIS lepton-proton 
scattering, the faster-than-linear growth of (n) with In 14^, a feature already well-known from 
e~^e~ annihilations and hadron-hadron collisions, and expected in perturbative QCD. 

Various models predict the evolution of the mean multiplicity with energy. We have fittedQ 
several parameterisations to the data plotted in Fig. 2a. 



According to the KNO-G prescription an appropriate power- law form is, following |T6| 

{n) = a- (W/Wo)^'' - c, (10) 

where the constant c may be regarded as a "discreteness correction" with the value c = 0.5 and 
Wo = l GeV. The fit yields a = 1.40 ± 0.04, b' = 0.20 ± 0.01 with xV^DF = 39/23. The value 
for b' agrees with that obtained by Szwed et al. {b' = 0.221) in a comprehensive analysis of the 
full phase space multiplicity distribution in e~^e~ annihilation . 



We have further compared the DIS data to the Modified Leading-Log (MMLA-(-LPHD) 
prediction in the form proposed in and valid for running QCD coupling 



{n) =ci^Nla + C2, (11) 



with 



NUY) = T{B){-) h^siz), (12) 



where z = JfY, with Y = \n{W/2Qo), a = 11 + 2A^//27, (3o = 11 - (2/3)A^/, B = a//3o 



and Nf the number of active flavours; I^, is a modified Bessel function of order u and F is the 
Gamma function. 



In this expression was shown to describe the mean charged multiplicity in e e~ annihi- 
lation from LEP energies down to centre-of-mass energies of 3 GeV. The factor | accounts for 
the multiplicity difference in a quark and gluon jet; Ci is a (non-perturbative) normalisation 
parameter and C2 the "leading parton" contribution, not included in the theoretical calcula- 



tion. Using the same shower cut-off value Qq = 270 MeV and Nf = 3 as in |£9[, we find 
ci = 1.21 ±0.05, C2 = 0.81 ±0.08 with xV^DF = 45/23. The best-fit curve is shown in Fig. 2a 
(dashed line) and describes the data over a wide W range. 

^For clarity, only a representative sample of vp data is plotted. 

^Because of the systematic difference between the EMC and E665 results, discussed in (29|, fits are performed 
without the E665 data points. The errors on best-fit parameters quoted in this section include systematic 
uncertainties. 
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The mean multiplicity has also been computed as a function of as{W) including the resum- 
mation of leading and next-to-leading corrections with the result 



{n) = aa^ exp (c/ y as) [1 + d ■ y/cQ , (13) 

where the parameter a cannot be calculated from QCD. The constants b and c are predicted 
by theory]]. This QCD prediction has been successfully tested over a wide energy range in 
several analyses of the mean charged multiplicity in e~^e~ annihilation including 



the recent LEP measurement at 130 GeV |Q. For the running coupling constant we use the 
two-loop expression 

as{W^) 1 /3ilnln(VrVA2) 



47r /3o\n{WyA^) f3l\n\wyA^) 



(14) 



with the constant Po as defined before, pi = 102 - 38A^//3, 6 = 1/4+ (10iV/)/(27/3o) and 
c = V967r/ Pq. According to [§0| A needs not be identical to Aj^, although both are expected 



to be rather similar, in particular if the 0(^/aJ) correction turns out to be small. 

A fit of the data to the form (^) with a and d as free parameters, A = 263 MeV and 



Nf = 3 yields a = 0.041 ± 0.006 and d = 0.2 ± 0.3 with xVnDF = 27/23. Neglecting the 
0{y/a^) correction in (|IB|) and treating A as a free parameter, we find A = 190 ± 60 MeV and 
a = 0.034 ±0.005 with xV^DF = 29.9/23. The functions (0), and (|13D are fairly similaif] 



and clear differences between them will become visible only at much larger energies. 

The above comparisons confirm and extend earlier indications from np interactions ^2 



that the rate of increase with energy of (n) in the current fragmentation region of DIS lepton- 
nucleon interactions is similar to that observed in e'^e~ annihilation in the presently covered 
energy range. 

The energy evolution of the mean multiplicity of partons emitted from a primary parton 
has been calculated for running as well as for fixed as [^|6^ . In the latter case the multiplicity 
rises as a power of the energy as in equation (0). For running coupling the growth is slower 
than any power of W, but faster than any power of In W, as in equations (|lT]) and ([T3|) . 

The various parameterisations discussed here show that, up to present energies, distinction 
between fixed and running as, based on measurements of (n), is still not possible. The data are 
sensitive, however, to soft gluon interference in QCD. Neglecting interference would increase 
the multiplicity anomalous dimension [Q, 7 = dln(n)/dlnPF^, by a factor ^/2. This has been 



shown ||6^ to be inconsistent with data for any reasonable value of A and confirmed for DIS 
in 1,0. 

In a strict sense the perturbative QCD results apply only to the soft component of the 
parton cascade emitted by a single quark or gluon jet. They are approximately valid for initial 
parton configurations such as those encountered in e~^e~ annihilation and in the simple quark- 
parton picture of DIS. However, the dynamics of DIS processes is more complex than in e~^e~ 
and differences of detail are to be expected |6^ . 



In Fig. 2a we compare the DIS multiplicity with single hemisphere results as expected from 
the JETSET generator in e~^e~ annihilation for light quark- ant iquark pairs (dotted line)g. The 
evolution with energy of (n) in e~^e~ is similar to that in DIS. The HI results are consistent 
with the presence of a small multiplicity excess in e'^e~ annihilation, relative to DIS, above 

^Note that formula ( p^ ) reduces to at large z. 

®The corresponding fitted curves practically coincide with the dashed line in Fig. 2a and are not shown for 
reasons of clarity. 

^The single hemisphere multiplicity distribution for e+e^ is calculated relative to the thrust axis. 
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W = 10 GeV, which has been noted before p5|,[55[]. In measurements of inclusive charged 
particle spectra at HERA [^,|68| a similar excess is seen near zero longitudinal momentum in 
the hadronic centre-of-mass. It is usually attributed to more prolific gluon emission in e~^e~. 
The MEPS 6.4 generator for DIS at the HERA energy (solid line) overestimates substantially 
the mean charged multiplicity. 



4.1.2 The ratio {n)/D and R2 

In Fig. 2b we show the ratio of (n) to in a comparison with fixed target DIS data. This ratio 
is expected to be energy independent if KNO scaling holds. The data above 10 GeV are indeed 
constant, within large errors, after a clear rise at lower W. The JETSET e~^e~ prediction 
(dotted curve) exhibits rather similar energy dependence and, moreover, illustrates that KNO 
scaling is only approximately valid. The MEPS 6.4 generator (solid line) agrees well with the 
HERA data. 

The normalised second order factorial moment R2 = {n{n — 1)) / (n)^ is plotted in Fig. 2c. 
This quantity is equal to the integrated two-particle inclusive density and is, therefore, a direct 
measure of the strength of hadron-hadron correlations. It shows little, if any, energy dependence 
over the HERA range but rises steadily at lower W. The behaviour is very similar for e~^e~ final 
states. The experimental values of R2 are further compared with a QCD calculation which, for 
a quark jet, predicts R2 to behave as 

7 

R2 = - [I — K^/a2 , (15) 
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with K = 0.88 for three flavours p4[. Leading and next-to-leading order predictions are plotted, 
with as calculated according to the two-loop formula (|T3D and the scale parameter set to 
A = 263 MeV. Both curves are significantly above the data. Nevertheless, it is interesting 
that the next-to-leading order calculation comes closer to the data, although the disagreement 
remains considerable. The data are rather well reproduced by the JETSET model in the case of 
e~^e~ annihilation, a fact confirmed by differential measurements of the two-particle correlation 



function by OPAL The prediction of MEPS 6.4 for DIS (solid line) is in agreement with 
the HERA results. 

The results presented in this and the previous section indicate similarities of the low order 
moments of the single hemisphere multiplicity distribution in DIS and e~^e~ annihilation, in 
conformity with the hypothesis of approximate universality of quark and gluon fragmentation. 

QCD predictions for the energy dependence of the mean parton multiplicity, derived from 
analytical solutions of the evolution equations, are in agreement with that observed for hadrons 
and thus add to the existing support for the LPHD ansatz at the single-inclusive level. The 



large disagreement between data and QCD next-to-leading order calculations of P3 for R2 



suggests that extension of the LPHD hypothesis to higher order inclusive correlations may not 
be justified [p3 . 



4.2 The multiplicity distribution in pseudorapidity domains 

In this section we study the multiplicity distribution and its moments in domains of pseu- 
dorapidity, limited to the current fragmentation region t]* > 1. We examine the and W 
dependence and compare with data from other types of interactions and with predictions from 
the MEPS 6.4 generator. 
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4.2.1 dependence 



In the simple quark-parton model the properties of the total hadronic system produced in a 
deep inelastic lepton-hadron collision depend on the lepton kinematical variables Bjorken-x 
and only through the invariant mass W of the hadronic system. In QCD, scaling violation 
of the quark fragmentation functions and of the parton distributions introduce an explicit 
dependence even at fixed W . Fixed target DIS electron, muon and (anti)neutrino experiments 
at low energies confirm that the global characteristics of the hadronic final states and the average 
number of produced hadrons in particular, vary most significantly with W . At fixed W , only 
weak dependences on are observed [0,0. Such results are in accord with the Bjorken-Kogut 
correspondence principle |7l| and imply that the densities in the hadron plateaus spanning the 
current, central and target regions are quite similar |]7T],|72|. Only recently has a statistically 
significant dependence of the mean charged hadron multiplicity been established in /x+p 
and v{v)p interactions JTSl, 173. The effect is limited to a restricted region in Feynman-x: 



—0.15 < xp < 0.15, where xp is the fractional longitudinal momentum of a hadron in the 
hadronic centre-of-mass frame. There are no published results on a possible variation with 
of the shape of the multiplicity distribution. 

The HI experiment allows us to investigate the charged particle multiplicity distribution over 
a widely extended range of in a novel energy domain. In Fig. 3 are plotted the mean charged 
multiplicity and the dispersion of the multiplicity distribution in four intervals of W , and in 
the pseudorapidity domain 1 < r/* < 5, covering part of the current fragmentation hemisphere. 
Within errors, no significant variation with Q^, in a fixed interval of W , is observed for (n) and 
D in the range 10 < < 1000 GeV^. The solid lines in Fig. 3 are fits to a constant. 

To ascertain the evolution with seen at much lower energy it will be necessary to study 
the multiplicity distribution in the central and proton remnant pseudorapidity regions which 
are not covered in this analysis. A similar study with quasi-real photons at HERA would be of 
evident interest. 

The data presented in subsequent sections have been obtained from data samples averaged 
over the full region (Q^ > 10 GeV^) covered by this experiment. 



4.2.2 The shape of the multiplicity distribution 

In Fig. 4 we show, as a representative example, the multiplicity distribution in the interval 
115 < W < 150 GeV, measured in various rj* domains. In this figure only statistical errors are 
plotted. The distribution for the widest rj* interval is plotted at its true scale. The distributions 
for the other intervals are successively shifted down by a factor ten. 

The figures illustrate that the multiplicity distribution, at fixed W , becomes narrower as 
the size of the rj* interval is reduced. However, the same distributions plotted in KNO form 
(not shown) widen under the same conditions. The latter property is in part related to the 
diminishing influence of global conservation constraints and was first predicted in 

Figure 4 also shows the predictions from MEPS 6.4 (open symbols). Significant deviations 
are noted. The model overestimates the mean multiplicity (cf. Fig. 2a) and, consequently, does 
not reproduce the small and high n tail of the distribution. These defects are seen in all W 
and Tj* intervals examined. 

Numerous parameterisations for the shape of the multiplicity distribution are used in the 
literature. Here we concentrate on the Negative Binomial Distribution (NDB) and the Lognor- 
mal Distribution (LND). The phenomenological arguments leading to these forms are discussed 
in Sect. 0. The parameters of these parametric models are obtained from a least fit to the 
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uncorrected multiplicity distributions, as explained in Sect. 3.4. The best-fit parameters for all 
studied pseudorapidity and W intervals are summarised in Table ||. The errors quoted are the 
quadratic sum of the statistical error and the systematic uncertainties. 

The solid (dashed) line in Fig. 4 shows how the LND (NBD) compares to the measurements. 
Inspection of this figure and of Table ^ indicates that the LND gives a reasonably accurate 
description of the data, in particular in the smallest 77* domain. However, the quality of the 
fits deteriorates in larger domains. Likewise, the NBD fits are acceptable in the smallest rj* 
domain but become progressively worse for larger intervals. The two distributions are seen 
to differ most for low multiplicities. Nevertheless, it may be verified from Tables 4 and 5 
that the estimates for the mean and dispersion of the multiplicity distribution, derived from 
these parameterisations, agree very well with those derived from the fully unfolded distribution. 
Parametric forms, such as the NBD and LND, therefore remain useful for phenomenology. 

In intermediate size rapidity intervals in e~^e~ annihilations at LEP both the LND and the 
NBD are unable to describe the multiplicity distribution which exhibits a prominent shoulder 
at intermediate n values. The shoulder is most prominent in the single hemisphere distri- 
butions [^, This structure results from a superposition of two-jet and three- or four-jet 
events. The effect demonstrates that the fluctuations in the number of hadrons, and therefore 
the multiplicity distribution, carries information on the hard partonic phase of the multihadron 
production process, even after soft hadronisation. The connection between parton level and 
hadron level dynamics has been investigated by the Lund group for parton cascades treated in 
the dipole formalism |]75| , [76[| . It is shown that, for centre-of-mass energies above ~ 50 GeV, the 
fluctuations in the hadron multiplicity are to better than 90% determined by the hardest and 
second hardest gluons emitted. Further softer radiation and subsequent (string) hadronisation, 
adds only small (sub-Poissonian) fluctuations to those induced in the initial stage of shower 
development. These analytical results provide a quantitative realisation of the notion of Local 
Parton-Hadron Duality derived directly from perturbative QCD. 

The HI measurements show no evidence for a shoulder structure of the type seen in e~^e~ . 
Neither is such a structure present in the multiplicity distribution predicted by MEPS 6.4. An 
excess of high-multiplicity events would also be expected if a significant proportion of DIS events 
were induced by QCD instantons In order to determine from our data an upper limit for 
the cross section of such events, we closely follow an analysis method recently applied by HI in a 



study of strange particle production ||5T|. We assume that the observed multiplicity distribution 
is a superposition of two distributions, one associated with instanton-induced events, another 
with "standard" DIS events. The former is calculated from the instanton generator described 
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5TI. For the latter we adopt a Negative Binomial form. Using a minimisation procedure 



to determine the relative proportion of instanton-induced and standard DIS events, we derive 
a 95% confidence level upper limit of 0.3 nb on the cross section for instanton production in 
the pseudorapidity domain 1 < 77* < 5 and 80 < < 115 GeV. 

The HI data, which cover larger W values than presently reachable in e~^e~ annihilation, are 
in qualitative accord with the QCD expectation that hard multi-jet production is less frequent 
in lepton-hadron collisions. Radiation from the target remnant — a composite colour source — 
is expected to be strongly suppressed by the so-called "antenna effect" |7BI]. For sufficiently 



small values of Bjorken-x it also suppresses hard gluon radiation in the current region. The 
probability per event for hard gluon radiation is consequently larger in e^e~ annihilation. As 
shown in [^], the antenna effect leads to smaller values of the local multiplicity anomalous 
dimensions and to smaller multiplicity fluctuations for DIS in the current region, compared to 
e~^e~ annihilation. These topics are examined in the following sections. 
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4.2.3 KNO scaling and correlations 



To demonstrate the energy scaling of the multiplicity distribution in DIS in the energy range 
opened up by HERA we show in Fig. 5 the KNO distributions "^{z) in the domain 1 < t]* < 5 
for four intervals in W, plotted on a logarithmic (top) and a linear scale (bottom), respectively. 
The dotted curve is the KNO function for e+e" annihilation[^ in the same i]* range, obtained 
from the JETSET parton shower model at a centre-of-mass energy of 91.2 GeV. The two data 
sets are remarkably similar. 

Exact KNO scaling implies that the function "^{z), and hence the moments Cq, are inde- 
pendent of W. Fig. 6 displays the variation of C moments with W for various t]* intervals (see 
also Table 4). At fixed W the moments increase as the rj* domain decreases in size, reflecting 
the widening of the multiplicity distribution in KNO form. 

The moments in the smallest and largest pseudorapidity domain show, within errors, little 
W dependence and thus exhibit approximate KNO scaling. However, violation of KNO scaling 
is seen in intermediate size intervals. This observation is consistent with the clear KNO scaling 
violations observed at HERA in DIS data on multiplicity distributions measured in the current 
region of the Breit frame of reference P,p!0[]. 

The MEPS 6.4 generator (solid line) describes the data well in the largest r/* domains but 
tends to underestimate C2 and C3 in the smaller ones. 

The dotted lines in Fig. 6 are expectations from JETSET for e~^e~ in the same pseudora- 
pidity domains as covered by HI. Large differences between DIS and e~^e~ annihilation are 
predicted, in particular for higher order moments in small 77* domains. However, in the largest 
domains the e^e~ results join smoothly with the DIS data. It will be of interest to confirm 
these predictions with future measurements at LEP in the W range studied here. 

Also shown is a measurement for non-single-diffractive pp collisions at ^/s = 200 GeVj^ in 
the interval 1 < |?7*| < 2 from UA5 [|1^. Here the fluctuations in particle density number near 



the central plateau are significantly larger than in DIS. It will, however, be shown in Sect. [4.2.4 
that the particle density itself is quite similar in the two processes. 

The values of the cumulants Kg in a given domain of phase space are a direct measure 
of the strength of "genuine" correlations among hadrons. Inspection of Table 4 shows that 
the three-particle correlation function is significantly different from zero only in the smallest 
interval 1 < r^* < 2, in accord with measurements for other types of interactions 

To study the two-particle correlation function more directly we present in Fig. 7 (n) and 
the second order factorial moment R2 = 1 + K2 in 77* domains. The mean charged multiplicity 
increases approximately as \nW for 1 < r^* < 2, i.e. near the central region, but faster in 
larger domains. The W dependence of R2 for 1 < r^* < 2 is less clear, in view of the errors, 
but compatible with the slow (logarithmic) rise well established in hadron-hadron and e+e" 
interactions. The MEPS 6.4 generator reproduces reasonably well the behaviour of R2 but 
systematically overestimates the mean multiplicity in all rj* domains. 



4.2.4 Particle density and Eq- flow 

The flow of transverse energy, Et, in multiparticle final states at high energy is studied in- 
tensively at hadron colliders and by the HERA experiments. The Et distribution in a given 

^°Single hemisphere e+e" data are known to exhibit KNO scaHng above w 20 GcV [|6|, |7) . The JETSET 
simulations show the same in the domain 1 < rj* < 5 discussed here. 

"The UA5 measurement at ^/s = 900 GeV (not shown) yields C2 = 1.84 ± 0.02, C3 = 4.0 ± 0.1 and 
C4 — 14.5 ± 0.8, the same, within errors, as the results at 200 GeV. 
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phase space domain, and its moments, are convoluted observables. They depend not only on 
the particle density and the pt structure of the collisions but also on multiparticle correlations, 
and, therefore, on the moments of the multiplicity distribution in that domain |7^. Here we 
compare the evolution with W of the particle density to that of the mean Et for DIS and 
hadron-hadron interactions. 

Fig. 8 shows a compilation of measurements in DIS P5| , ^ and in non-single-diffractive 



hadron-hadron collisions |T^,^|8^,^,q^ of mean charged multiplicity per unit of rapidity or 

pseudorapidity (solid symbols). The HI data are the same as those shown in Fig. 7a. The 

open symbols show measurements of the mean Et in the region —0.5 <•)]*< 0.5, presented 

and discussed in |^3|. The solid line shows a parameterisation used by UAl |^ with the form 
(n) = 0.35 + 0.74 (iy2)0.i05_ 

In spite of the differences in the rapidity region covered, the known difference between 
pseudorapidity and rapidity density, different experimental procedures and systematics, it re- 
mains of interest to note that the charged particle density for DIS at HERA interpolates quite 
smoothly with DIS and hadron-hadron data at lower and much higher energy. This observation 
is consistent with the analogy between virtual photon-hadron, real photon-hadron and hadron- 
hadron interactions originally advocated by Gribov and Feynman[^ which suggests universality 
of dynamics in the central plateau [^. However, the difference in correlation strength, noted 
in the previous section, suggests that such "universality", while applicable to single-inclusive 



spectra [£3[, may not hold for higher order correlations. 

Several mechanisms have been suggested to explain the dynamics of the large multiplicity 
fluctuations in soft hadron-hadron collisions: impact parameter averaged Poisson-like fluctu- 
ations; multiple soft parton interactions in the same event leading to mini-jets [Q; multi- 



pomeron exchange as in the Dual Parton Model |^5[ and multi-string conflgurations as in the 



Lund FRITIOF model [g6|. These mechanisms have no direct analogues in DIS, except for 
the boson-gluon fusion QCD process which can lead to two-string colour topologies and could, 
therefore, mimic multi-string properties in hadron-hadron models. 

In Fig. 8 we further compare the energy evolution of (n) in the central region in DIS with 
that in e~^e~ annihilation. The dotted curve is the prediction from the JETSET generator for 
the interval 1 < rj* < 2. It shows that the hadron density evolves much faster with W than in 
the DIS and hadron-hadron data. 

The energy evolution of (n) in perturbative QCD is controlled by the anomalous multi- 
plicity dimension 7 through 7 = dln(n)/dlniy^. For e~^e~ annihilation, and in restricted 



(pseudo)rapidity intervals, it is given by 7ee('7*) ~ y 3as(fcTmax)/27r, where kTmax is the max- 
imum possible transverse momentum at a given rj* [|76|. From the JETSET e^e^ predictions 
above W = 20 GeV we derive that 7ee is constant with a value of 0.16. This is somewhat smaller 



than expected from the analytic result (~ 0.22 at 200 GeV) but agrees with that quoted in 

In leptoproduction, where gluon emission from the nucleon remnant is thought to be sup- 
pressed, 'jDisiv*) depends more strongly on rj* with 'jDisiv*) ~ (V2)7ee('7*) for not too large 
positive 7]* From the DIS data in Fig. 8, using the UAl parameterisation, we estimate 

7z)/s = 0.8 — 0.9, consistent with expectations. This is the flrst semi- quantitative experimental 
conflrmation of the "antenna suppression" effect in DIS. However, better data are needed, also 
at lower W, to exploit these perturbative QCD predictions in a fully quantitative manner. This 
result also implies that the rate of increase with energy of (n) in the full current hemisphere 
for DIS and that in a single hemisphere for e+e~ annihilations, discussed in Sect. 4.1.1, are 

-^^For further discussion and references see fSSl], Chapt. 4. 
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expected to differ at higher centre of mass energies. 

Non-asymptotic analytical QCD predictions for higher order multiplicity moments are at 
present not available. It is known, however, that the energy dependence of the moments is 
asymptotically controlled by •yoisiv*) leeiv*) iHl- It is therefore likely that the differences 



between the C moments in DIS and e^e~, described in Sect. [4.2.3| , are a further reflection of 
suppressed gluon emission in DIS. 

Finally, the comparison between (n) and {E^) in Fig. 8 demonstrates that both evolve 
with energy rather similarly. This suggests that the increase of mean transverse energy with 



increasing W (decreasing xb) at fixed Q , previously observed in this experiment follows 



mainly from an increase of the hadron multiplicity and less so from a rise of the mean transverse 
momentum of individual hadrons. A direct measurement of the energy dependence of the 
transverse momentum distribution in or near the central region of deep inelastic ep collisions at 
HERA, and comparison with existing hadron-hadron data, should help to clarify this interesting 
and theoretically much debated issue. 



5 Summary 

Data, fully corrected for detector effects, are presented on the evolution with W and of 
the charged particle multiplicity distribution and its statistical moments, over the ranges 80 < 
W < 220 GeV and 10 < < 1000 GeV^ in subdomains of pseudorapidity space, including 
the full current hemisphere. The main results can be summarised as follows: 

• The mean charged hadron multiplicity and the dispersion, measured in fixed intervals of 
W and in the domain 1 < r^* < 5 show, within errors, no dependence on the virtuality of 
the exchanged boson over the range covered by HI. 

• The low order moments of the multiplicity distribution in the full current hemisphere show 
noteworthy similarities with single hemisphere data in e~^e~ annihilation, in conformity 
with the hypothesis of approximate environmental independence of quark hadronisation. 
In particular, the mean charged hadron multiplicity in DIS shows a similar rate of increase 
with W to that measured in e~^e~ annihilation up to 130 GeV. Analytical predictions 
from perturbative QCD on the mean parton multiplicity in jets, which are proven to be 
valid for hadrons in e+e~ annihilation, are therefore confirmed for the first time in DIS 
at higher energies than presently reachable at LEP. Data on the second order factorial 
moment show that higher order QCD corrections and/or non-perturbative effects remain 
significant at present energies. 

• The analysis of the multiplicity distribution in pseudorapidity domains of varying size 
proves that the well-documented property of KNO scaling, a general characteristic of 
scale-invariant stochastic branching processes, remains valid in DIS at HERA for small 
and for large pseudorapidity intervals, but not in intermediate size domains. The KNO 
phenomenon, also predicted in QCD at asymptotic energies, results from an intricate 
interplay of correlations, different in different types of collisions, and changing rapidly 
over phase space. The KNO function in the region 1 < r^* < 5 is strikingly similar to that 
expected for e+e~ annihilation under the same kinematical conditions. 

• The charged particle density near the central region of the 7*p centre-of-mass system 
grows significantly more slowly in DIS than in e^e~ annihilation. The strength of particle 
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correlations, as reflected in the C moments, is much larger in the latter process for small 
pseudorapidity domains. Such differences can be understood within perturbative QCD 
from calculations of the local anomalous multiplicity dimensions within the Lund dipole 
formalism. The comparison between e^e~ annihilation and DIS data provides direct 
evidence for the "antenna suppression" effect in current fragmentation, a characteristic 
of deep inelastic lepton-hadron dynamics at small Bjorken-a;. The same mechanism offers 
a qualitative explanation for the absence of a multi-jet induced shoulder structure in 
the DIS multiplicity distributions in intermediate size pseudorapidity domains, now well 
estabhshed in e+e~ annihilation data at LEP. 

• The charged particle density near the central plateau in the deep inelastic process is 
of the same magnitude as that for minimum bias non-single-diffractive hadron-hadron 
interactions at the same value of the centre-of-mass energy. Its evolution with W is 
also comparable to that in hadron-hadron collisions measured up to 900 GeV. However, 
hadron-hadron collisions are characterised by substantially stronger correlations in small 
pseudorapidity domains. These features remain to be understood within the Gribov- 
Feynman pictures of DIS. 

• A comparison of the evolution with W of mean transverse energy flow and charged particle 
density near the central region in DIS shows that the two phenomena are strongly cor- 
related. The striking similarity with the behaviour seen in non-single-diffractive hadron- 
hadron collisions implies that a theoretical explanation within QCD should simultaneously 
address the dynamics of both types of process. 

• The multiplicity distributions in the smallest r)* domain examined can be well parame- 
tcriscd with Lognormal or Negative Binomial functions. However, the quality of the fits 
deteriorates in larger domains. The large n tail of the experimental distributions is well 
described but deviations occur for small multiplicities. Nevertheless, as economic repre- 
sentations of the data, these parametric forms continue to be useful for phenomenology. 

• Among a variety of Monte Carlo generators presently being developed for DIS, we have 
shown predictions from MEPS 6.4. This model overestimates the mean charged multi- 
plicity in the full current hemisphere, as well as in subdomains of pseudorapidity space. 
Higher order moments of the multiplicity distributions are reasonably well described. 
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W range (W) range (Q^) % events 
(GeV) (GeV) (GeV^) (GeV^) 
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201.9 
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13.9 
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40 


27.6 
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54.9 


1037 
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374.2 


349 








10 -> 


80 


23.5 


9255 



Table 1: Number of events, mean Q'^ and mean W in the different kinematical regions studied. 



W (GeV) 


80 115 


r/* > 

115 150 150 4- 185 


185 -4- 220 


80 + 115 


1< ry* < 3 

115 + 150 150 + 185 


185 + 220 


photo- 


-0.7% 


-0.8% 


-0.8% 


-0.0% 


-0.9% 


-1.2% 


-0.5% 


+0.3% 


production 


(1.3%) 


(1.5%) 


(1.6%) 


(2.2%) 


(1.7%) 


(2.4%) 


(2.5%) 


(3.3%) 


BEMC ener- 


-2.1% 


-1.7% 


-0.8% 


-hO.2% 


-0.4% 


-1.1% 


+2.7% 


+2.4% 


gy scale 


(1.2%) 


(1.5%) 


(1.5%) 


(1.7%) 


(1.4%) 


(2.0%) 


(2.2%) 


(2.8%) 


track 


-4.0% 


-3.7% 


-3.9% 


-3.1% 


-3.0% 


-3.5% 


-4.3% 


-4.5% 


selection 


(1.2%) 


(1.5%) 


(1.6%) 


(1.5%) 


(2.1%) 


(1.9%) 


(2.0%) 


(2.4%) 


unfolding 


-^0.1% 


-0.8% 


+0.9% 


+1.8% 


+0.6% 


-1.2% 


+1.4% 


+1.1% 


method 


(1.2%) 


(1.4%) 


(1.7%) 


(1.8%) 


(1.7%) 


(1.9%) 


(2.4%) 


(2.7%) 


LEPTO 6.1 


-4.9% 


-4.5% 


-4.6% 


-3.8% 


-0.7% 


-2.3% 


-2.5% 


-2.1% 




(1.6%) 


(0.9%) 


(1.2%) 


(1.2%) 


(1.1%) 


(1.0%) 


(1.2%) 


(1.2%) 


TOTAL 


6.8% 


6.6% 


6.6% 


5.9% 


3.7% 


5.3% 


6.4% 


6.5% 



Table 2: Summary of systematic effects. The change in average multiphcity is given for the full 
current hemisphere and for a restricted pseudorapidity interval. For comparison, the statistical 
error on (n) for the given data sample is given in brackets. 
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80 GeV < < 115 GeV 



n P„ (%) 

1< r?* < 2 1< ??* < 3 1< r?* < 4 1< t?* < 5 






15.79 


± 


0, 


,52 


± 


2.18 


3.26 


± 


0.22 


± 


1, 


,58 


0.78 


± 


0.09 


± 


1.66 


0, 


.40+ ± 


0.08 








1 


22.55 


± 


0, 


,50 


± 


1.92 


8.01 


± 


0.27 


± 


1, 


,29 


2.73 


± 


0.18 




0.94 


1, 


.64 


± 


0.16 


± 


0, 


,39 


2 


20.62 


± 


0, 


,46 


± 


1.24 


12.84 


± 


0.35 


± 


1, 


,30 


5.73 


± 


0.28 


± 


1.45 


3, 


.82 


± 


0.22 


± 


1, 


,27 


3 


15.96 




0, 


,43 




2.29 


14.03 




0.35 


± 


1, 


,21 


9.55 


± 


0.29 


± 


1.46 


7, 


,59 


± 


0.28 


± 


1, 


,06 


4 


10.21 




0, 


,37 


± 


1.06 


13.38 




0.32 


± 


2, 


,22 


11.11 


± 


0.32 




1.21 


10, 


,08 


± 


0.34 


± 


1, 


,18 


5 


6.07 


± 


0, 


,30 


± 


1.54 


11.73 


± 


0.35 


± 


0, 


,46 


12.76 


± 


0.35 


± 


2.22 


12, 


.30 


± 


0.34 


± 


0, 


,89 


6 


3.85 


± 


0, 


,26 


± 


0.48 


9.84 


± 


0.34 


± 


0, 


,71 


12.19 


± 


0.32 


± 


1.28 


12, 


.91 


± 


0.35 


± 


1, 


,29 


7 


2.24 


± 


0. 


,16 


± 


0.62 


7.94 


± 


0.33 


± 


0, 


,52 


11.17 


± 


0.37 


± 


0.96 


12, 


.27 


± 


0.35 


± 


1, 


.86 


8 


1.15 


± 


0. 


.12 


± 


0.42 


5.91 


± 


0.32 


± 


1, 


,20 


9.22 


± 


0.31 


± 


0.48 


10 


.65 


± 


0.33 


± 


0, 


.86 


9 


0.68 


± 


0. 


,10 


± 


0.31 


4.00 


± 


0.21 


± 


0, 


,82 


7.72 


± 


0.33 


± 


1.37 


8, 


.69 


± 


0.30 


± 


1, 


.08 


10 


0.49 


± 


0. 


,12 


± 


0.24 


3.21 


± 


0.19 


± 


0, 


,56 


5.69 


± 


0.24 


± 


1.31 


6, 


.38 


± 


0.26 


± 


0, 


.99 


11 


0.19 


± 


0. 


,06 




0.26 


2.22 


± 


0.18 


± 


1, 


,06 


3.95 


± 


0.22 


± 


1.00 


4, 


.79 


± 


0.25 


± 


1, 


.23 


12 


0.07 


± 


0. 


,04 




0.13 


1.49 




0.14 


± 


0, 


,66 


2.76 


± 


0.23 


± 


0.96 


3, 


.20 


± 


0.24 


± 


0, 


.97 


13 


0.05 


± 


0. 


,03 


± 


0.06 


0.90 


± 


0.15 


± 


0.53 


1.79 


± 


0.18 


± 


0.80 


2, 


.05 


± 


0.19 


± 


0, 


.69 


14 














0.54 


± 


0.11 


± 


0, 


,53 


1.22 


± 


0.15 


± 


0.62 


1, 


.38 


± 


0.16 


± 


0, 


,74 


15 














0.36 




0.08 


± 


0, 


,57 


0.67 


± 


0.12 


± 


0.68 


0, 


,86 


± 


0.14 


± 


0, 


,66 


16 














0.14 


± 


0.07 


± 


0, 


,19 


0.40 


± 


0.11 


± 


0.43 


0, 


,40 


± 


0.08 


± 


0, 


,24 


17 














0.08 


± 


0.05 


± 


0, 


,10 


0.24 


± 


0.11 


± 


0.33 


0, 


.27 


± 


0.13 


± 


0, 


,37 


18 


























0.13 


± 


0.08 


± 


0.19 


0, 


.10 


± 


0.05 


± 


0, 


,11 



115 GeV <W <150 GeV 

n Pn (%) 

Kij* <2 1< T?* < 3 1< ?;* < 4 1< ?;* < 5 






16. 


,80 


± 


0.62 


± 


0.55 


4.52 


± 


0.35 


± 


1. 


1 


21. 


,55 


± 


0.54 


± 


1.05 


8.00 


± 


0.35 


± 


1. 


2 


19. 


,69 




0.50 


± 


0.53 


12.15 




0.36 


± 


0. 


3 


15, 


,19 


± 


0.48 


± 


0.77 


12.96 


± 


0.36 


± 


0, 


4 


10, 


,54 


± 


0.36 


± 


0.31 


12.74 


± 


0..38 


± 


0, 


5 


6, 


,76 


± 


0.34 


± 


0.75 


11.09 


± 


0.32 


± 


1, 


6 


4, 


17 


± 


0.27 


± 


0.50 


9.35 


± 


0.31 


± 


0, 


7 


2, 


,55 


± 


0.20 


± 


0.76 


7.31 


± 


0.27 


± 


2, 


8 


1, 


,29 


± 


0.14 


± 


0.48 


6.01 


± 


0.29 


± 


0, 


9 


0, 


,64 


± 


0.10 


± 


0.62 


4.53 




0.24 


± 


1, 


10 


0, 


,31 


± 


0.07 


± 


0.14 


3.76 


± 


0.25 


± 


0, 


11 


0, 


,20 




0.04 




0.07 


2.57 


± 


0.23 


± 


0, 


12 


0, 


12 




0.03 


± 


0.15 


1.90 


± 


0.18 


± 


0, 


13 


0, 


,07 


± 


0.03 


± 


0.07 


1.05 


± 


0.16 


± 


0, 


14 


0, 


,03 


± 


0.02 


± 


0.04 


0.77 


± 


0.12 


± 


0, 


15 














0.56 


± 


0.10 


± 


0, 


16 














0.31 


± 


0.08 


± 


0, 


17 














0.18 


± 


0.06 


± 


0, 


18 














0.12 


± 


0.06 


± 


0, 



19 
20 
21 



19 


1.09 


± 


0.18 


± 


0.77 


0.35+ ± 


0.08 






21 


2.57 


± 


0.19 


± 


1.12 


1.17 


± 


0.13 


± 


0.62 


70 


4.72 


± 


0.24 


± 


0.70 


2.58 


± 


0.14 


± 


1.26 


91 


8.09 


± 


0.29 


± 


0.74 


5.63 


± 


0.23 


± 


0.69 


74 


9.96 


± 


0.33 


± 


0.95 


7.97 




0.24 


± 


1.03 


63 


11.23 


± 


0.36 


± 


0.95 


10.36 


± 


0.31 


± 


0.82 


70 


11.27 


± 


0.33 


± 


1.37 


11.44 


± 


0.37 


± 


0.56 


06 


10.32 


± 


0.31 




0.25 


11.96 


± 


0.35 


± 


1.18 


93 


9.21 


± 


0.34 


± 


0.52 


10.41 


± 


0.32 


± 


0.66 


38 


8.03 


± 


0.30 


± 


0.29 


9.65 


± 


0.36 


± 


0.74 


53 


6.51 




0.29 


± 


0.87 


7.88 




0.28 


± 


0.49 


86 


5.32 


± 


0.26 


± 


0.88 


6.15 




0.23 


± 


0.80 


77 


3.84 


± 


0.23 


± 


0.88 


4.96 


± 


0.27 


± 


1.35 


29 


2.96 


± 


0.23 


± 


0.53 


3.51 


± 


0.25 


± 


0.50 


54 


1.81 




0.18 




0.65 


2.15 




0.21 


± 


0.45 


48 


1.25 


± 


0.16 


± 


0.33 


1.47 


± 


0.17 


± 


0.43 


24 


0.76 


± 


0.12 


± 


0.62 


0.93 


± 


0.14 


± 


0.92 


28 


0.46 


± 


0.10 


± 


0.29 


0.72 


± 


0.17 


± 


0.62 


12 


0.25 


± 


0.09 


± 


0.20 


0.25 


± 


0.08 


± 


0.34 




0.14 


± 


0.06 


± 


0.23 


0.17 


± 


0.07 


± 


0.24 




0.08 


± 


0.05 


± 


0.11 


0.12 


± 


0.08 


± 


0.13 














0.03 


± 


0.02 


± 


0.03 



Table 3: The fully corrected multiplicity distribution P„ (%). The first error is the statistical 
error, the second the systematic uncertainty of the result. 

The value of Pq in the domain 1 < < 5 is not measured but taken from the reweighted 
DJANGO 6.0 Monte Carlo generator. 
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150 GeV <W <185 GeV 



n Pn (%) 

1< ?7* < 2 1< ?7* < 3 1< ?7* < 4 1< ??* < 5 






16.01 ± 0.58 ± 5.43 


4.24 


± 


0.36 


± 


3. 


,19 


1.13 


± 


0. 


,15 


± 


1, 


,51 


0.28t ± 


0. 


,08 






1 


19.56 ± 0.63 ± 1.10 


8.32 


± 


0.41 


± 


1. 


.39 


2.56 


± 


0. 


,20 


± 


1, 


,55 


1.09 


± 


0. 


,13 


± 


1.03 


2 


19.73 ± 0.54 ± 1.58 


10.99 


± 


0.40 


± 


0. 


,67 


4.69 


± 


0. 


,24 


± 


2, 


.38 


2.21 


± 


0, 


,14 


± 


0.71 


3 


16.01 ± 0.60 ± 1.37 


11.75 


± 


0.37 


± 


1. 


,30 


6.81 


± 


0, 


,30 


± 


1, 


.00 


4.18 


± 


0, 


,20 


± 


1.28 


4 


11.19 ± 0.50 ± 1.20 


11.71 


± 


0.41 


± 


1. 


,27 


8.77 


± 


0, 


,35 


± 


1, 


.02 


6.61 


± 


0, 


,23 


± 


2.08 


5 


7.21 ± 0.34 ± 0.70 


11.37 


± 


0.41 


± 


1. 


,49 


9.87 


± 


0, 


,38 


± 


1, 


.56 


8.24 


± 


0, 


,31 


± 


1.36 


6 


4.11 ± 0.30 ± 0.81 


9.39 


± 


0.39 


± 


1, 


,08 


10.47 


± 


0, 


,34 


± 


0, 


,32 


9.77 


± 


0, 


,34 


± 


1.20 


7 


2.80 ± 0.25 ± 0.93 


7.98 


± 


0.35 


± 


0, 


,57 


9.68 


± 


0, 


,33 


± 


0, 


,94 


11.04 


± 


0, 


,33 


± 


0.89 


8 


1.40 ± 0.17 ± 0.81 


6.27 


± 


0.38 


± 


0, 


,52 


9.10 




0, 


,34 


± 


0, 


,57 


10.87 


± 


0, 


,34 


± 


0.66 


9 


0.91 ± 0.15 ± 0.48 


5.36 


± 


0.35 


± 


1, 


,07 


8.10 


± 


0, 


,34 


± 


1, 


,05 


9.89 


± 


0, 


,35 


± 


1.13 


10 


0.33 ± 0.08 ± 0.23 


3.64 


± 


0.25 


± 


1, 


,03 


6.90 


± 


0, 


,31 


± 


0, 


,81 


8.18 


± 


0, 


,35 


± 


0.44 


11 


0.31 ± 0.07 ± 0.18 


2.87 


± 


0.22 


± 


0, 


,28 


5.76 


± 


0, 


,30 




0, 


,72 


7.07 


± 


0, 


,31 




0.56 


12 


0.10 ± 0.04 ± 0.10 


1.82 


± 


0.20 




0, 


,73 


4.92 


± 


0, 


,28 




1, 


,17 


5.76 




0, 


,34 




1.15 


13 


0.13 ± 0.07 ± 0.18 


1.24 


± 


0.15 


± 


0, 


,54 


3.59 


± 


0, 


,25 




1, 


,32 


4.88 


± 


0, 


,31 


± 


1.50 


14 


0.08 ± 0.06 ± 0.14 


1.08 


± 


0.16 


± 


0, 


,89 


2.71 


± 


0, 


,24 


± 


0, 


,97 


3.42 


± 


0, 


,26 


± 


0.89 


15 




0.78 


± 


0.15 


± 


0, 


,61 


1.59 


± 


0, 


,16 


± 


0, 


,72 


2.43 


± 


0, 


,25 


± 


1.00 


16 




0.43 


± 


0.09 


± 


0, 


,40 


1.11 


± 


0, 


,12 


± 


0, 


,93 


1.31 


± 


0, 


,17 


± 


0.78 


17 




0.20 


± 


0.06 


± 


0. 


,13 


0.86 


± 


0, 


,13 


± 


0.46 


1.07 


± 


0, 


,19 


± 


0.64 


18 




0.18 


± 


0.07 


± 


0. 


,13 


0.49 


± 


0, 


,13 


± 


0, 


.32 


0.55 


± 


0, 


,14 




0.51 


19 




0.13 


± 


0.07 


± 


0. 


,24 


0.31 


± 


0, 


,11 


± 


0, 


.63 


0.43 


± 


0, 


,13 


± 


0.50 


20 




0.06 


± 


0.05 


± 


0. 


,10 


0.09 


± 


0, 


,05 


± 


0, 


.14 


0.12 


± 


0, 


,05 


± 


0.51 


21 
















0.17 


± 


0, 


,09 


± 


0, 


.30 


0.18 


± 


0, 


,08 


± 


0.32 


22 
















0.13 


± 


0, 


,09 


± 


0, 


.26 


0.16 


± 


0, 


,11 


± 


0.31 



185 GeV <W < 220 GeV 



1< r?* < 2 



Pn (%) 

1 < < 3 1 < ??* < 4 



1 < r;* < 5 






15.76 


± 


0.68 




3.11 


4.24 




0.42 


± 


3.18 


1 


19.66 


± 


0.76 


± 


0.78 


8.00 


± 


0.45 


± 


0.48 


2 


19.56 


± 


0.62 


± 


1.75 


10.74 


± 


0.57 


± 


1.52 


3 


15.76 


± 


0.61 


± 


2.03 


11.51 


± 


0.47 


± 


0.69 


4 


11.08 


± 


0.58 


± 


1.85 


11.76 


± 


0.48 


± 


2.08 


5 


6.95 




0.45 


± 


0.86 


11.65 




0.53 


± 


1.27 


6 


4.63 


± 


0.38 


± 


0.97 


9.74 


± 


0.47 


± 


1.62 


7 


3.20 


± 


0.32 


± 


1.14 


7.96 


± 


0.43 


± 


0.91 


8 


1.36 


± 


0.16 


± 


0.74 


6.83 




0.44 


± 


1.06 


9 


0.95 


± 


0.15 




0.79 


5.12 




0.34 


± 


1.16 


10 


0.47 




0.12 


± 


0.31 


3.70 


± 


0.30 


± 


1.24 


11 


0.33 


± 


0.10 


± 


0.50 


2.52 


± 


0.24 


± 


0.88 


12 


0.11 




0.04 


± 


0.09 


1.79 


± 


0.19 




0.76 


13 


0.05 


± 


0.04 


± 


0.10 


1.39 


it 


0.17 


± 


0.92 


14 


0.03 


± 


0.02 


± 


0.05 


1.11 


± 


0.16 


± 


1.23 


15 












0.66 


± 


0.13 


± 


0.33 


16 












0.44 




0.10 


± 


0.27 


17 












0.29 


± 


0.10 


± 


0.19 


18 












0.14 


± 


0.06 


± 


0.19 


19 












0.10 


± 


0.05 


± 


0.14 


20 












0.06 


± 


0.04 


± 


0.29 


21 












0.05 


± 


0.04 


± 


0.10 



22 



1.07 


± 


0.24 


± 


2, 


,96 


0.27t ± 


0.09 






2.81 


± 


0.30 


± 


1, 


,27 


0.96 


± 


0.13 


± 


1.96 


5.22 


± 


0.32 


± 


0, 


,73 


2.37 


± 


0.23 


± 


0.68 


6.70 


± 


0.33 


± 


0, 


,50 


3.73 


± 


0.24 


± 


1.12 


7.98 


± 


0.37 


± 


0, 


.87 


5.47 


± 


0.27 


± 


0.99 


9.18 


± 


0.34 


± 


0, 


.58 


7.50 


± 


0.30 


± 


0.83 


9.07 


± 


0.40 


± 


0, 


.87 


8.49 


± 


0.35 


± 


1.24 


10.18 


± 


0.50 


± 


0, 


.63 


10.08 


± 


0.39 


± 


1.61 


9.10 


± 


0.41 


± 


0, 


.68 


10.54 


± 


0.46 


± 


1.18 


7.95 


± 


0.34 


± 


1, 


.57 


10.40 


± 


0.43 


± 


1.71 


7.50 


± 


0.34 


± 


1, 


.83 


9.20 


± 


0.44 


± 


1.00 


6.10 


± 


0.38 


± 


0.78 


8.03 


± 


0.40 


± 


1.12 


4.98 




0.33 


± 


1, 


,30 


6.48 


± 


0.30 


± 


0.99 


3.46 


± 


0.26 


± 


0, 


,77 


4.78 


± 


0.29 


± 


0.95 


2.60 


± 


0.22 


± 


0, 


,75 


3.81 


± 


0.27 


± 


1.22 


2.03 


± 


0.21 


± 


0, 


,58 


2.32 


± 


0.21 


± 


0.78 


1.29 




0.16 




0, 


,80 


2.13 




0.22 


± 


1.08 


1.00 


± 


0.16 


± 


0, 


,65 


1.14 


± 


0.14 


± 


0.67 


0.71 


± 


0.15 


± 


0, 


,22 


0.83 


± 


0.12 


± 


0.44 


0.39 


± 


0.09 


± 


0, 


,35 


0.64 


± 


0.14 


± 


0.29 


0.21 


± 


0.07 


± 


0, 


,18 


0.30 


± 


0.07 


± 


0.39 


0.17 


± 


0.10 


± 


0, 


.26 


0.13 


± 


0.07 


± 


0.18 


0.11 


± 


0.04 


± 


0, 


.13 


0.19 


± 


0.08 


± 


0.26 














0.03 


± 


0.02 


± 


0.05 



Table 3: continued 
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r)* >0 

W (GeV) 80 115 115 4- 150 150 185 185 4- 220 



{rich} 


8.98 


± 0.07 


± 0.61 


10.00 


± 


0.07 


± 0.66 


10, 


.88 


± 


0.09 


± 0.72 


11.35 


± 


0.09 


± 


0.67 


D2 


3.89 


± 0.07 


± 0.23 


4.16 


± 


0.06 


± 0.26 


4, 


,49 


± 


0.08 


± 0.36 


4.56 


± 


0.08 


± 


0.30 


Da 


3.43 


± 0.16 


± 0.22 


3.41 


± 


0.15 


± 0.21 


3, 


,61 


± 


0.26 


± 0.55 


3.67 


± 


0.25 


± 


0.45 


D4 


5.39 


± 0.15 


± 0.38 


5.58 


± 


0.13 


± 0.23 


6, 


,05 


± 


0.21 


± 0.71 


6.26 




0.19 


± 


0.49 


Ca 


1.18 


± 0.00 


± 0.04 


1.17 


± 


0.00 


± 0.01 


1, 


,17 


± 


0.00 


± 0.02 


1.16 


± 


0.00 


± 


0.02 


Cs 


1.62 


± 0.02 


± 0.14 


1.56 




0.01 


± 0.06 


1, 


,54 


± 


0.02 


± 0.08 


1.51 


± 


0.02 


± 


0.07 


C4 


2.48 


± 0.07 


± 0.40 


2.29 


± 


0.05 


± 0.17 


2, 


,26 


± 


0.07 


± 0.21 


2.19 


± 


0.06 


± 


0.19 


R2 


1.07 


± 0.00 


± 0.03 


1.07 


± 


0.00 


± 0.01 


1, 


,07 


± 


0.00 


± 0.02 


1.07 


± 


0.00 


± 


0.02 


Rs 


1.24 


± 0.02 


± 0.10 


1.22 


± 


0.01 


± 0.04 


1, 


,24 


± 


0.02 


± 0.07 


1.22 


± 


0.02 


± 


0.06 



Krj* <2 

W (GeV) 80 -r 115 115 -4- 150 150 -4 185 185 4- 220 



{rich} 


2.46 


± 0, 


,03 


± 0.10 


2.50 


± 


0.03 


± 


0.12 


2, 


,63 


± 


0.04 


± 0.18 


2.66 


± 


0.04 


± 


0. 


,18 


D2 


2.07 


± 0, 


,05 


± 0.08 


2.13 


± 


0.03 


± 


0.10 


2, 


,22 


± 


0.05 


± 0.12 


2.22 


± 


0.03 


± 


0. 


14 


Da 


2.21 


± 0, 


11 


± 0.10 


2.24 


± 


0.05 


± 


0.14 


2, 


,39 


lb 


0.11 


± 0.18 


2.31 


± 


0.07 


± 


0, 


,12 


Di 


3.08 


± 0, 


16 


± 0.15 


3.15 


± 


0.07 


± 


0.17 


3, 


,39 


± 


0.16 


± 0.26 


3.24 


± 


0.09 




0, 


,12 


C2 


1.71 


± 0, 


,02 


± 0.02 


1.72 


± 


0.02 


± 


0.04 


1, 


,71 


± 


0.02 


± 0.05 


1.69 


± 


0.02 


± 


0, 


,03 


Ch 


3.87 


± 0, 


,16 


± 0.18 


3.90 




0.11 


± 


0.24 


3, 


,89 


± 


0.17 


± 0.18 


3.75 


± 


0.12 


± 


0, 


,20 


Ca 


10.67 


± 0, 


,93 


± 1.23 


10.77 


± 


0.60 


± 


1.17 


11, 


,03 


± 


1.09 


± 0.63 


10.02 




0.65 




1, 


,22 


R2 


1.30 


± 0, 


,02 


± 0.03 


1.32 


± 


0.01 


± 


0.04 


1, 


,33 


± 


0.02 


± 0.04 


1.32 




0.02 


± 


0, 


,03 


Ra 


2.11 


± 0, 


,14 


± 0.15 


2.15 


± 


0.08 


± 


0.19 


2, 


,22 


± 


0.14 


± 0.14 


2.11 


± 


0.09 


± 


0. 


13 


Ks 


0.187 


± 0, 


,070 ± 0.097 


0.167 ± 


0.040 ± 


0.063 


0, 


,226 


± 


0.079 ± 0.085 


0.147 ± 


0.043 ± 


0. 


,101 



1< r?* < 3 

W (GeV) 80 ^ 115 115 -r 150 150 4 185 185 -r 220 



{rich} 


4.90 


± 0. 


,04 


± 0.18 


5.06 


± 


0.05 


± 


0.27 


5, 


,32 


± 


0.07 


± 0.34 


5.35 


± 


0.06 


± 


0.35 


D2 


3.10 


± 0, 


,04 


± 0.13 


3.37 


± 


0.04 


± 


0.21 


3, 


,57 


± 


0.06 


± 0.26 


3.58 


± 


0.07 


± 


0.30 


D,. 


2.91 


± 0, 


,10 


± 0.25 


3.18 


± 


0.09 


± 


0.21 


3, 


,45 


± 


0.12 


± 0.43 


3.54 


± 


0.20 




0.34 


Di 


4.27 


± 0, 


,11 


± 0.28 


4.62 


± 


0.10 


± 


0.27 


5, 


,01 


± 


0.13 


± 0.45 


5.18 


± 


0.28 




0.45 


C2 


1.40 


± 0, 


,01 


± 0.03 


1.44 


± 


0.01 


± 


0.02 


1, 


,45 


± 


0.01 


± 0.02 


1.44 


± 


0.01 


± 


0.02 


Ch 


2.41 


± 0, 


,04 


± 0.16 


2.57 


± 


0.05 




0.10 


2, 


,62 


± 


0.07 


± 0.14 


2.62 


± 


0.09 


± 


0.09 


Ca 


4.83 


± 0, 


,19 


± 0.64 


5.34 


± 


0.20 


± 


0.38 


5, 


,57 


± 


0.28 


± 0.68 


5.71 


± 


0.47 


± 


0.34 


R2 


1.19 


± 0, 


,01 


± 0.03 


1.24 


± 


0.01 


± 


0.02 


1, 


,26 


± 


0.01 


± 0.02 


1.26 


± 


0.01 


± 


0.03 


Ra 


1.64 


± 0, 


,04 


± 0.13 


1.80 


± 


0.04 


± 


0.09 


1, 


,87 


± 


0.06 


± 0.13 


1.88 


± 


0.08 


± 


0.11 


Ka 


0.048 


± 0, 


,015 ± 0.044 


0.063 


± 


0.016 ± 


0.033 


0, 


,088 


± 


0.022 ± 0.077 


0.108 


± 


0.041 ± 


0.036 



1 < ?7* < 4 

W (GeV) 80 -4- 115 115 ~ 150 150 4- 185 185 4- 220 



{rich} 


6.45 


± 0.06 


± 0.33 


7.00 


± 


0.05 


± 0.35 


7.51 


± 0.06 


± 


0.51 


7.66 


± 


0.07 


± 0.47 


D2 


3.23 


± 0.09 


± 0.14 


3.57 


± 


0.04 


± 0.23 


3.93 


± 0.06 


± 


0.25 


4.06 


± 


0.06 


± 0.20 


Da 


2.71 


± 0.28 


± 0.21 


2.89 


± 


0.11 


± 0.23 


3.28 


± 0.14 


± 


0.50 


3.32 


± 


0.16 


± 0.28 


Di 


4.34 


± 0.27 


± 0.24 


4.72 


± 


0.10 


± 0.20 


5.29 


± 0.14 


± 


0.44 


5.41 


± 


0.14 


± 0.22 


C2 


1.25 


± 0.01 


± 0.03 


1.26 


± 


0.00 


± 0.02 


1.27 


± 0.00 


± 


0.02 


1.28 


± 


0.00 


± 0.02 


Ca 


1.82 


± 0.05 


± 0.13 


1.85 




0.02 


± 0.07 


1.90 


± 0.03 


± 


0.08 


1.92 




0.03 


± 0.08 


C4 


3.01 


± 0.20 


± 0.38 


3.04 




0.08 


± 0.22 


3.22 


± 0.11 


± 


0.27 


3.26 


± 


0.12 


± 0.26 


i?2 


1.09 


± 0.01 


± 0.03 


1.11 


± 


0.00 


± 0.02 


1.14 


± 0.00 


± 


0.02 


1.15 


± 


0.00 


± 0.02 


Ra 


1.29 


± 0.05 


± 0.09 


1.35 


± 


0.02 


± 0.06 


1.43 


± 0.03 


± 


0.07 


1.45 


± 


0.03 


± 0.06 



Table 4: Moments and cumulants of the unfolded multiplicity distribution for different ran] 
in pseudorapidity rf and intervals in W . 
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1 < r/* < 5 

W (GeV) 80-^ 115 115 ^ 150 150 ~ 185 185 ^ 220 





6.90 


± 0.05 


± 0.34 


7.72 


± 


0.05 


± 0.42 


8 


.45 


± 0, 


,06 


± 0.58 


8.81 


± 


0.08 


± 


0.54 


D'z 


3.15 


± 0.05 


± 0.21 


3.45 


± 


0.04 


± 0.25 


3 


.77 


it 0, 


,06 


± 0.25 


3.86 


± 


0.05 


± 


0.22 


D-A 


2.56 


± 0.12 


± 0.33 


2.74 


± 


0.13 


± 0.19 


3 


.02 


± 0, 


,18 


± 0.54 


2.95 


± 


0.19 


± 


0.41 


Di 


4.23 


± 0.11 


± 0.37 


4.59 


± 


0.10 


± 0.29 


5 


.08 


± 0, 


,15 


± 0.45 


5.17 


± 


0.14 


± 


0.29 


Ci 


1.20 


± 0.00 


± 0.03 


1.19 


± 


0.00 


± 0.02 


1 


.19 


± 0. 


.00 


± 0.02 


1.19 


± 


0.00 


± 


0.02 


Cs 


1.67 


± 0.02 


± 0.12 


1.64 


± 


0.02 


± 0.06 


1 


.64 


± 0. 


.02 


± 0.08 


1.61 


± 


0.02 


± 


0.08 


C4 


2.60 


± 0.07 


± 0.35 


2.50 


± 


0.06 


± 0.16 


2 


.50 


± 0. 


.08 


± 0.24 


2.42 


± 


0.07 


± 


0.21 


Ri 


1.06 


± 0.00 


± 0.02 


1.07 


± 


0.00 


± 0.02 


1 


.08 


± 0. 


.00 


± 0.02 


1.07 


± 


0.00 


± 


0.02 


Rz 


1.19 


± 0.02 


± 0.09 


1.21 


± 


0.01 


± 0.05 


1 


.24 


± 0. 


.02 


± 0.07 


1.23 


± 


0.02 


± 


0.06 



2 < r/* < 3 

W (GcV) 80 H- 115 115 H- 150 150 H- 185 185 ^ 220 





2.46 


± 


0.03 


± 0.17 


2.61 


± 


0.03 


± 


0.20 


2, 


,73 




0, 


,04 


± 0.27 


2.71 


± 


0.04 




0.20 


r>2 


1.92 


± 


0.03 


± 0.19 


2.11 




0.03 




0.16 


2, 


,24 


± 


0, 


,04 


± 0.20 


2.22 


± 


0.03 


± 


0.18 




1.86 




0.05 


± 0.10 


2.07 


± 


0.06 


± 


0.14 


2, 


,23 




0, 


,08 


± 0.13 


2.28 




0.07 


± 


0.22 


D4 


2.67 




0.06 


± 0.12 


2.94 


± 


0.08 


± 


0.21 


3, 


,14 


± 


0, 


,09 


± 0.18 


3.20 


± 


0.09 


± 


0.25 


C2 


1.61 


± 


0.01 


± 0.03 


1.65 


± 


0.01 


± 


0.03 


1, 


,67 


± 


0. 


,01 


± 0.05 


1.67 


± 


0.02 


± 


0.06 


Cz 


3.26 




0.08 


± 0.16 


3.47 


± 


0.10 


± 


0.15 


3. 


.57 


± 


0. 


,10 


± 0.25 


3.60 


± 


0.11 


± 


0.30 


Ci 


7.77 




0.35 


± 0.79 


8.57 


± 


0.48 


± 


0.60 


8, 


.99 


± 


0. 


.52 


d= 1.20 


9.30 


± 


0.56 


± 


1.19 


R2 


1.20 




0.01 


± 0.02 


1.27 


± 


0.01 


± 


0.04 


1. 


.30 


± 


0. 


,01 


± 0.05 


1.30 


± 


0.02 


± 


0.06 


Rs 


1.63 


± 


0.05 


± 0.10 


1.86 


± 


0.08 


± 


0.13 


1, 


.99 


± 


0. 


,08 


± 0.14 


2.02 


± 


0.09 


± 


0.25 


Ka 


0.018 


± 


0.021 


± 0.041 


0.043 ± 


0.030 


± 


0.035 


0, 


.072 


± 


0. 


.040 ± 0.062 


0.121 


± 


0.038 


± 


0.076 



3 < r?* < 4 

W (GeV) 80 ^ 115 115 -r 150 150 4- 185 185 4- 220 



{rich) 


1.54 


± 0, 


,02 


± 0.18 


1.93 


± 


0.02 


± 


0.17 


2, 


,17 


± 


0.02 


± 0.18 


2.29 




0.03 




0.17 


D2 


1.32 


± 0, 


,01 


± 0.14 


1.48 


± 


0.02 


± 


0.14 


1, 


,59 


± 


0.02 


± 0.09 


1.67 


± 


0.03 




0.10 


Da 


1.23 


± 0, 


,03 


± 0.07 


1.33 


± 


0.03 


± 


0.09 


1, 


,39 


± 


0.05 


± 0.04 


1.54 


± 


0.05 


± 


0.16 


D4 


1.79 


± 0, 


,03 


± 0.11 


2.01 


± 


0.04 


± 


0.10 


2, 


,14 


± 


0.06 


± 0.06 


2.31 


± 


0.06 


± 


0.17 


C-z 


1.73 


± 0, 


,01 


± 0.11 


1.58 


± 


0.01 




0.10 


1, 


,53 


± 


0.01 


± 0.06 


1.53 


± 


0.01 


± 


0.10 


Ca 


3.70 


± 0. 


,08 


± 0.55 


3.08 


± 


0.07 


± 


0.47 


2, 


.87 


± 


0.07 


± 0.26 


2.91 


± 


0.08 


± 


0.46 


Ca 


9.24 


± 0. 


.38 


± 2.29 


6.99 


± 


0.28 


± 


1.72 


6, 


.21 


± 


0.31 


± 1.00 


6.47 


± 


0.34 


± 


1.70 


R2 


1.08 


± 0. 


,01 


± 0.04 


1.07 


± 


0.01 


± 


0.07 


1, 


.07 


± 


0.01 


± 0.03 


1.09 


± 


0.01 


± 


0.08 


Rs 


1.17 


± 0. 


,05 


± 0.15 


1.16 


± 


0.04 


± 


0.19 


1, 


.18 


± 


0.05 


± 0.08 


1.28 


± 


0.05 


± 


0.24 


Ka 


0.079 


± 0. 


,019 ± 0.062 


0.045 ± 


0.014 ± 


0.040 


0, 


.054 


± 


0.018 ± 0.033 


0.011 


± 


0.017 ± 


0.068 



4 < r/* < 5 

W (GcV) 80 H- 115 115 ^ 150 150 H- 185 185 4- 220 



(rich) 


0.46 


± 


0.01 


± 


0.04 


0.71 




0.01 




0.04 


0.93 


± 


0.01 


± 


0.09 


1.17 


± 


0.01 




0.09 


D2 


0.71 


± 


0.01 


± 


0.02 


0.88 


± 


0.01 


± 


0.08 


0.99 


± 


0.01 


± 


0.06 


1.09 


± 


0.01 


± 


0.19 


Da 


0.83 




0.02 


± 


0.02 


0.93 


± 


0.01 




0.03 


0.98 




0.02 


± 


0.05 


1.01 


± 


0.02 


± 


0.01 


Di 


1.09 


± 


0.03 


± 


0.04 


1.26 


± 


0.02 


± 


0.03 


1.37 




0.03 


± 


0.07 


1.47 


± 


0.02 


± 


0.03 


C-z 


3.34 


± 


0.06 


± 


0.11 


2.52 


± 


0.04 


± 


0.18 


2.13 


± 


0.02 


± 


0.09 


1.86 


± 


0.01 


± 


0.11 


Ca 


13.78 


± 


0.61 


± 


1.06 


7.85 


± 


0.30 


± 


1.18 


5.58 


± 


0.17 


± 


0.59 


4.23 


± 


0.10 


± 


0.55 


Ci 


68.80 


± 


5.70 


± 


9.55 


28.93 


± 


1.85 


± 


6.83 


17.21 


± 


0.87 


± 


3.47 


11.24 


± 


0.46 


± 


2.25 


R2 


1.19 


± 


0.04 


± 


0.06 


1.12 


± 


0.03 




0.10 


1.05 


± 


0.02 


d= 


0.05 


1.01 


± 


0.01 


± 


0.05 


Ra 


1.41 


± 


0.24 


± 


0.29 


1.16 


± 


0.09 


± 


0.23 


1.02 


± 


0.07 


± 


0.19 


0.92 


± 


0.04 


± 


0.09 


Ka 


0.159 ± 


0.144 ± 


0.138 


0.221 ± 


0.041 ± 


0.136 


0.150 ± 


0.028 ± 


0.104 


0.109 ± 


0.017 ± 


0.061 



Table 4: continued 
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1< r/* < 2 



W (GeV) 


80 115 


115 150 


150 185 


185 -f 220 


n 


2.52 ± 0.10 


2.58 ± 0.13 


2.72 it 0.15 


2.75 ± 0.19 


1/k 


0.285 ± 0.080 


0.288 ± 0.056 


0.269 it 0.034 


0.270 ± 0.064 


X / ' J— ' -T 


29.4/13 


23.0/14 


'^'^ 0/1 ^ 


27.1/14 


771 


9 -t- n 1 


9 4Q 4- n 1 R 


9 f\(\ -t- n 90 
z.uo in u.^u 


^.ui HI U.J-O 


d 


0.791 ± 0.092 


0.812 ± 0.077 


0.743 ± 0.094 


0.768 it 0.122 




2.6 ± 3.4 


3.1 ± 1.2 


1.9 it 1.7 


2.5 it 7.5 


X^/NDF 


21.8/12 


8.9/13 


24.9/14 


19.8/13 






1 < ??* < 3 






W TGeVI 

V V I V_J \3 V } 


80 -r 115 


115 150 


150 185 


185 4- 220 


n 


4.85 ± 0.16 


5.12 ± 0.26 


5.36 lb 0.25 


5.41 it 0.20 


1/k 


0.161 ± 0.039 


0.235 it 0.017 


0.256 lb 0.028 


0.245 it 0.034 




191 R /I Q 


o^.o/ ly 




44 7 /24 




4 7q -1- n 94 


o.Kjo HZ u.ou 




r: QO -1-0 17 
0.0^ in u. J. f 


d 


0.634 ± 0.017 


0.644 ± 0.080 


0.662 it 0.074 


0.641 it 0.060 




4.5 ± 2.2 


2.0 it 6.5 


2.1 it 5.5 


1.8 it 6.3 


X^/NDF 


71.1/18 


42.7/18 


43.6/20 


31.1/23 






1 < r/* < 4 






W (GeV\ 

V V \ V_X v> V ) 


80 115 


115 ~ 150 


1 50 — 1 85 


185 H- 220 


n 


6.41 ± 0.28 


7.00 it 0.39 


7.52 zb 0.46 


7.71 it 0.30 


1/k 


0.085 ± 0.018 


0.102 ± 0.014 


0.132 lb 0.026 


0.142 it 0.020 


v2 /AJ n F 


Qfi S /91 


112.4/21 


oo.u / zo 


44.7/24 






u. yo iL u. o / 




7 c;q J- n oq 
1 .Oi? ~ u. ^0 




0.475 ± 0.009 


0.484 ± 0.011 


0.505 it 0.023 


0.527 ± 0.022 


c 


2.7 ± 1.3 


2.7 it 2.5 


2.8 it 7.1 


4.5 it 5.9 


X^/NDF 


62.6/20 


69.6/20 


43.8/24 


32.8/23 






1 < ??* < 5 






W (GeV) 


80 -r 115 


115 -=-150 


150 ~ 185 


185 -f- 220 


n 


6.90 ± 0.33 


7.73 ± 0.41 


8.44 ± 0.49 


8.88 ± 0.34 


l/fe 


0.067 ± 0.011 


0.068 it 0.012 


0.074 it 0.013 


0.069 it 0.020 


X^/NDF 


44.7/19 


57.0/21 


80.8/24 


63.6/25 


m 


6.89 ± 0.29 


7.72 it 0.41 


8.40 it 0.48 


8.78 it 0.33 


d 


0.433 ± 0.010 


0.422 it 0.007 


0.422 it 0.015 


0.422 ± 0.020 


c 


1.7 ± 1.5 


1.8 ± 2.0 


2.4 it 5.1 


4.1 ± 6.1 


X^/NDF 


41.6/18 


52.8/20 


58.4/23 


27.7/24 



Table 5: Parameters of Negative Binomial and Lognormal distributions in domains of pseudo- 
rapidity T)* and W. Errors quoted are the quadratic sum of the statistical error and systematic 
uncertainties. 



32 



° 2 




A 
C 

V 



a 



10 



• H1 e> 

■ E665 |iD 
^ EMCjip 
▼ WA25vp 

• WA21 vp 

• JETSETeV 
- MEPS6.4 




MLLAfit ^ w twj^ 



A 
C 

V 



b) 



2.5 
2 
1.5 



! I I 



1 + 



T i 


i 


i ! 


• 











1.75 -• 



1.5 



c. 



1.25 



leoding order 



next— to-leading order^_ _ 



0.75 



10 



10 



W (GeV) 



Figure 2: (a) W dependence of the mean charged particle multiphcity in the full current 
hemisphere in comparison with fixed target lepton-nucleon data. The various curves are de- 
scribed in the text; (b) and (c) same as (a) for the ratio {n) / D and the second factorial 
moment R2, respectively. 
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Figure 3: The Q'^ dependence of the mean charged particle multiphcity (top) and the dispersion 
(bottom) in intervals of W and in the domain 1 < 77* < 5. Horizontal lines are fits to a constant. 
The inner error bars represent the statistical error, the outer error bars represent the total 
(quadratic sum of statistical and systematical) errors. 
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Figure 4: The unfolded multiplicity distribution in the interval 115 GeV < W < 150 GeV 
and > 10 GeV^, in indicated pseudorapidity domains. The distribution for 1 < r^* < 5 
is plotted at its true scale; each consecutive distribution is shifted down by factor of 10. The 
HI data points (solid symbols) are compared to MEPS 6.4 (open symbols) and to fits with 
the Lognormal (full) and Negative Binomial (dashed) distribution. Statistical errors only are 
plotted. 
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Figure 5: Unfolded multiplicity distributions in the pseudorapidity domain 1 < 77* < 5 mea- 
sured in different W intervals plotted in KNO form, in logarithmic and linear scale. Statistical 
errors only are plotted. 
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Figure 6: W dependence of the multiplicity moments Cg in various rj* domains. Data from 
UA5 1 19 1 in the interval 1 < \ri*\ < 2 are also shown. From top to bottom, the domains are 
1 < ?7* < 2, 1 < 77* < 3, 1 < ?7* < 4 and 1 < rj* < 5. The curves are described in the text. 
The inner error bars represent the statistical errors, the outer error bars represent the total 
(quadratic sum of statistical and systematical) errors. 
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Figure 7: W dependence of mean charged multiplicity and second order factorial moment R2 
in indicated pseudorapidity domains, compared with MEPS 6.4 predictions. The inner error 
bars represent the statistical errors, the outer error bars represent the total (quadratic sum of 
statistical and systematical) errors. 
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Figure 8: Solid symbols: W dependence of charged particle density in DIS and hadron- 
hadron collisions (left scale); open symbols: the mean transverse momentum flow, per unit of 
(pseudo) rapidity (right scale). The curves are described in the text. 77* is pseudorapidity, while 
y* is rapidity in the hadronic centre-of-mass frame. 
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